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ABSTRACT 

We present a CARMA 1.3 mm continuum survey toward 9 Class 0 protostars in the 
Perseus molecular cloud at ~0 / . / 3 (70 AU) resolution. This study approximately doubles 
the number of Class 0 protostars observed with spatial resolutions < 100 AU at millime¬ 
ter wavelengths, enabling the presence of protostellar disks and proto-binary systems to 
be probed. We detect flattened structures with radii > 100 AU around 2 sources (L1448 
IRS2 and Per-emb-14) and these sources may be strong disk candidates. Marginally- 
resolved structures with position angles within 30° of perpendicular to the outflow are 
found toward 3 protostars (L1448 IRS3C, IRAS 03282+3035, and L1448C) and are 
considered disk candidates. Two others (L1448 IRS3B and IRAS 03292+3039) have 
resolved structure, possibly indicative of massive inner envelopes or disks; L1448 IRS3B 
also has a companion separated by Cf/9 (~210 AU). IC348-MMS does not have well- 
resolved structure and the candidate first hydrostatic core L1451-MMS is marginally 
resolved on 1" scales. The strong disk candidate sources were followed-up with C 18 0 
(J = 2 —> 1) observations, detecting velocity gradients consistent with rotation, but it 
is unclear if the rotation is Keplerian. We compare the observed visibility amplitudes 
to radiative transfer models, finding that visibility amplitude ratios suggest a compact 
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component (possibly a disk) is necessary for 5 of 9 Class 0 sources; envelopes alone may 
explain the other 4 systems. We conclude that there is evidence for the formation of 
large disks in the Class 0 phase with a range of radii and masses dependent upon their 
initial formation conditions. 

Subject headings: planetary systems: proto-planetary disks — stars: formation — tech¬ 
niques: interferometric — stars: protostars 


1. Introduction 


Stars form due to the gravitational collapse of dense condensations within molecular clouds and 
the star formation process is what ultimately leads to the formation of solar systems and potentially 
life itself, as has transpired in our Solar System. Therefore, the formation of the proto-planetary 
disk is a key step in the process of forming a planetary system. Moreover, disks may also play a 


role in the formation of binary and multiple star systems via gravitational instability (e.g., Kratter 


et al. 2010); about one-third of all star systems and 50% of sun-like star systems are found in binary 
or multiple systems with a typical separation of ~50 AU (Raghavan et al. 2010). 


The formation of disks is thought to begin early in the star formation process, during the Class 


0 phase of protostellar evolution (Andre et al. 1993). This is the earliest recognizable phase of the 


star formation process, characterized by a protostar surrounded by a dense envelope of gas and 
dust. The early formation of disks is inferred from the ubiquity of outflows and jets during the 
protostellar phase (e.g. Frank et al.|2 014, and references therein). Disks form during the collapse of 
the star forming cloud via conservation of angular momentum (e.g. [Ulrich 1976; Cassen &; Moosman 
1981, Terebey et al. 1984); the angular momentum may derive from bulk cloud rotation or the net 
angular momentum of residual turbulent motions. 

Disk formation, however, may not be as simple if the protostellar cloud is magnetized (even 
weakly) ( Allen et al.||2003 Mellon Sz Li|2008 Li et al.||2013 ). The magnetic field will be dragged in¬ 
ward with the collapsing material, increasing the local field strength. The magnetic field then slows 
the rotation of the inner envelope by transporting angular momentum to the larger-scale envelope 


(Allen et al. 2003. Mellon Sz Li 2008). This ‘magnetic-braking’ can slow the rotation efficiently 


enough to suppress disk formation and this result has become known as the ‘magnetic-braking 
catastrophe.’ More recently, a number of ways around the catastrophe have been found: reduction 


of magnetic held strengths via Ohmic dissipation (Dapp Sz Basu 2010), misaligned rotation axes 


and magnetic held directions ( Joos et aL||2012 ), turbulence (Seifried et al. 2013), reconnection dif¬ 
fusion (Lazarian et al.|[2012), and/or reduced ionization (hence less held coupling) (Padovani et al. 


2013). 


While the theoretical difficulty of disk formation is resolved somehow by nature, Class 0 disks 
remain poorly-characterized observationally. There have been efforts to detect and resolve Class 0 
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disks with millimeter interferometers for the past 20 years. Possible signatures of disks have been 


detected via dust continuum emission (e.g., 

Chandler et al. 

1995; 

Brown et al. 2000 

Looney et al. 

2000. 

Harvey et al. 

2003, 

Jprgensen et al. 

2009 

); however, these pioneering efforts were not able 


to uniquely identify disks due to lack of sensitivity and/or spatial resolution. Further muddying 


the waters, Maury et al. (2010) did not resolve structure consistent with disks having R > 100 
AU, using data from the Plateau de Bure Interferometer (PdBI) with ~0 / . / 5 resolution. However, 


the Maury et al. (2010) sample was comprised of two Very Low Luminosity Objects (VeLLOs) 
(IRAM 04191 and L1521F, Dunham et al. 2006 Bourke et al. 2006) and three typical Class 


0 sources (L1527 IRS, L1448C, NGC 1333 IRAS2A). Another high-resolution study (0 / . / 3) with 
the Combined Array for Research in Millimeter-wave Astronomy (CARMA) toward the Class 0 


protostellar system L1157-mm also failed to resolve a disk (Chiang et al. 2012). 


Following these studies, the Class 0 protostar L1527 IRS (also in the Maury et al. sample) was 


found to have evidence for an edge-on disk in 3.8 /mi scattered light (Tobin et al. 2010a). Subsequent 
follow-up observations of this source led to the disk being resolved in the dust continuum as well 


as confirmation of rotational support (Tobin et al. 2012, 2013). The mass of the protostar was 
measured to be ~0.2 M & , about 20% of the surrounding envelope mass. Observations of L1527 
with the Atacama Large Millimeter/submillimeter Array (ALMA) now show that the mass of the 
protostar is likely 0.3 Mq, with a rotationally-supported radius of 54 AU (Ohashi et al.|2 01 1). Sakai 
et al. (2014bja) also found evidence for a chemical change at the interface between the envelope 
and disk using new data from ALMA. Evidence of rotationally-supported Class 0 disks in other 
systems have now been found with ALMA (Murillo et al.||20l3 Lindberg et al.|[2014 Codella et al. 


2014). Thus, the combined results of previous studies point toward a diversity of disk properties 


in the Class 0 phase and that large disks may form in some systems before the end of the Class 0 
phase. 

Despite these previous efforts, the number of Class 0 protostars observed with spatial resolu¬ 
tions better than 100 AU is only ~10 at millimeter and submillimeter wavelengths. To further the 
characterization of Class 0 disks, we have observed a sample of 9 Class 0 protostars and 2 Class I 
protostars in the Perseus molecular cloud (D ~ 230 pc; Hirota et al. 2008, 2011) using the CARMA 


array at 1.3 mm with a resolution of ~0 / . , 3 (70 AU). The Perseus molecular cloud was chosen for 
this study because it is the nearest star forming region with a large number of Class 0 protostars; 
the more nearby Taurus and Ophiuchus clouds do not have sufficiently large numbers of Class 0 
protostars from which to draw a meaningful sample. 


In this paper, we present the results from our survey, examining the dust continuum structures, 
visibility amplitudes, and compare to radiative transfer models. We will also present molecular line 
follow-up obtained toward some sources using CARMA and the Submillimeter Array (SMA) as well 
as the CO (J = 2 —> 1) outflow emission when detected. We present the sample, observations, and 
data reduction in Section 2, the continuum and molecular line results are described in Section 3, 
and a comparison to radiative transfer models is presented in Section 4. The results are discussed 
in Section 5 and we present our summary and conclusions in Section 6. 
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Observations 


2.1. Sample Selection 

We have drawn our sample of nine Class 0 sources (Table 1) from the known protostars in 
the Perseus molecular clouds based on the various infrared and millimeter studies that have been 
conducted (e.g., Hirano et al.||1999 Looney et al.||2000 Enoch et al.||2009 Chen et aT7||2010 Pineda 


et al. 2011, Schnee et al. 2012 Sadavoy et al. 2014). In addition, two Class I protostars were 


located adjacent to Class 0 sources and within the CARMA field of view. Enoch et al. (2009) 


classifies 27 sources as Class 0 protostars; however, additional sources are now known from more 


recent millimeter observations (references above) and one source in the Enoch et al. (2009) survey 


was incorrectly associated with IRAC emission (Per-emb-37). Including, these additional sources 
brings the total number of Class 0 protostars in Perseus to ~38. 

We chose sources that were not known to be close binaries (< 1000 AU) and had clearly defined 


outflow axes from CO observations or Spitzer IRAC imaging (e.g. Jprgensen et al. 2006). We also 
focused on sources that had not been previously observed with sensitive, high-resolution observa¬ 
tions, aiming to increase the number of sources with high resolution observations. Good knowledge 
of the outflow axis is important because this is a proxy for the current angular momentum vector 
of the system. If a sufficiently large disk is present, it is expected to be elongated perpendicular to 
the outflow direction. 


The Class 0 sources range in luminosity from 0.7 L 0 to 9.2 L 0 (Enoch et al. 2009, Sadavoy 
et al.|2014 ) and the candidate FHSC L1451-MMS has a luminosity less than 0.05 L 0 (Pineda et al 


2011). This range of luminosity is representative of the distribution of protostellar luminosities in 


Perseus, with only a few systems having luminosities in excess of 10 L 0 . The sources also sample 
a range of environments: wide binary systems, isolated cores within the cloud, members of the 
NGC1333 cluster, the outskirts of the older IC348 cluster, and the moderately active star forming 
cloud L1448. This sample enables us to begin more thoroughly exploring the presence of disks from 
their dust continuum emission. 


2.2. CARMA Observations and Data Reduction 

CARMA is a heterogeneous interferometer array located in the Inyo mountains of California. 
It is comprised of 6 x 10.4 m antennas, 9 x 6.1 m antennas, making up the main array, and 8 
x 3.5 m antennas that are operated as an independent array or as part of the CARMA23 array. 
The observations were conducted in B and C configurations of the main array between December 
2012 and November 2013, see Table 2. B configuration has baseline lengths up to ~1 km and 
C configuration has baseline lengths up to ~ 0.35 km. The respective angular resolutions of the 
configurations are ~0 / /3 and ~0 / . / 9. Three science targets were observed in each track, including a 
test source 0326+277 to verify the millimeter seeing in B configuration. 
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The observations were conducted in a standard loop, bracketing the three science pointings 
(and test source in B-configuration) with an observation of the phase calibrator (0237+288). Flux 
(Uranus) and bandpass (3C84) calibrators were observed at the beginning of each track. The 
C-PACS system had been used in previous extended configurations for correction of rapid phase 
variations (Perez et al. 2010), but the system was not available in the 2012 and 2013 observing 
seasons. However, the data were taken during exceptionally good weather conditions and standard 
calibration methods produced excellent data. The optical depth at 1.3 mm was typically ~ 0.1 
during B configuration and the root-mean-squared (RMS) path length (phase noise) was less than 
100 /irn as measured by the phase monitor (see Table 2); 100 /im RMS path length corresponds to 
~0 / . / 34 seeing, approximately the angular resolution provided by B-configuration at 1.3 mm. 


The central frequency of our observations was 225.0491 GHz (1.3 mm) and the correlator 
was configured for single-polarization observations with six 500 MHz sub-bands set for continuum 
observation, yielding 6 GHz of continuum bandwidth (dual-side band). The remaining two spectral 
windows were configured for spectral line observation with 31 MHz bandwidth and were set to 
observe the (J = 2 —> 1) transitions of 12 CO, 13 CO, and C 18 0; the individual channels were 97 
kHz corresponding to ~0.13 km s -1 velocity resolution. Follow-up observations of Per-emb-14 in 
November 2013 used dual-polarization mode in order to have higher spectral line sensitivity. 


We reduced and edited the visibility data using standard methods within the MIRIAD software 
package (Sault et al. 1995); see Tobin et al. (2013) or Chiang et al. (2012) for more details. The 
absolute flux calibration error is estimated to be ~20%. The data were imaged using the CLEAN 
algorithm using robust weighting to optimize the angular resolution and sensitivity. In order to 
ensure that our data reduction results in reliable images, we imaged the test source 0326+277 for 
each B-configuration data set. Gaussian fits to the source demonstrate that it is consistent with a 
point source, verifying the good millimeter seeing. 


2.3. SMA Observations and Data Reduction 

The SMA is an 8-element interferometer array comprised of 6.1 m antennas and is located on 
Mauna Kea in Hawaii. We observed the source L1448 IRS2 from the CARMA sample with the 
SMA in the Extended configuration on two dates in September 2013 (Table 3). We observed in 
two receiver mode with the low frequency receiver tuned to 225.434 GHz (1.3 mm) and the high 
frequency receiver tuned to 351.135 GHz (850 /im). In this mode, the correlator provides 4 GHz 
of bandwidth (2 GHz in each of the upper and lower sidebands) from each receiver. Each 2 GHz 
band is comprised of 24 sub-bands with a bandwidth of 104 MHz and our setup had 64 channels 
per sub-band by default and we used high-resolution settings to observe C 18 0 (J = 2 —> 1) (1024 
channels), 13 CO (J = 2 -+ 1) (512 channels) and 12 CO (J = 2 -+ 1) (512 channels) with the 
low frequency receiver. With the high frequency receiver, we observed the H 13 CO + and HCO + 
(J = 4 —> 3) transitions with 512 channels. 
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The visibility data were edited and calibrated using the MIR software package, an IDL-based 
software package originally developed for the Owens Valley Radio Observatory and adapted by the 
SMA group. The error in the absolute flux calibration is estimated to be ~10%. The calibrated 
visibilities were exported to MIRIAD format for imaging. We detected all the CO (J = 2 —> 1) 
isotopologues and HCO + (J = 4 —> 3), but we will only discuss the CO isotopologues and continuum 
results in this paper. 


2.4. Spatial Scale Sensitivity 

The incomplete sampling of the uv-plane limits the maximum size of structures that can be 
fully recovered with an interferometer. Moreover, the inclusion of only B and C or Extended and 
Compact configurations of CARMA and the SMA, respectively, limits the sensitivity to spatially 
extended objects. In order to characterize the sensitivity to different spatial scales, we generated 
a series of uniform brightness model images of symmetric disks having a variety of diameters and 
sampled them with the uv-coverage of the same configurations as our observations. We find that 
50% of the model flux density can be recovered at 1.3 mm from a 4" diameter disk in CARMA C- 
configuration, a E/75 disk for CARMA B-configuration, a 9 /7 disk for SMA Compact configuration, 
and a 4" disk for SMA Extended configuration. The SMA Extended configuration observations at 
850 fim can recover a 2f.'7 disk. Note that these sizes are optimistic because the uv-coverage in the 
simulations is not split between sources as in the observations and some observed data are flagged. 
Furthermore, a uniform brightness disk is not representative of the structure of protostellar sources, 
but this example is meant to present a limiting case. 


3. Results 

3.1. Continuum Emission Structure 

We detect 1.3 mm dust continuum emission toward all observed sources in both B and C con¬ 
figurations. The detected continuum structures broadly fit into four categories: disk-like structures 
(Section 3.1.1), candidate multiples (Section 3.1.2), featureless or unresolved (Section 3.1.3), and 
asymmetric resolved structures (Section 3.1.4). Some of the continuum data are highly structured 
in the combined B and C-array imaging, and others may only have resolved structure in the highest 
resolution data (B-array) alone. Images of the Class 0 protostars with sensitivity to structure on 
arcsecond scales with the combined B and C-array datasets are shown in Figure [lj The highest 
resolution images, emphasizing structure on sub-arcsecond scales, are shown in Figure [2] Further¬ 
more, two Class I protostars were within the primary beams of Class 0 sources and their images are 


°Per-emb-24 was only observed in C-array and was not observed in B-array due to only having a ~3<r (3.2 
mjy/beam) detection. 
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shown in Figure [3j The 1.3 mm flux densities, RMS noise, and resolution for each image are given 
Table 4. In addition, we fit two-dimensional Gaussians to the images shown in Figures 1 and 2. 
This was done with the imfit task in CASA to determine the deconvolved source sizes and position 
angles, these values are listed in Table 4. 


For a source to be classified as ‘disk-like,’ it must have a deconvolved position angle within 30° 
of orthogonal to the outflow axis; the relative position angles of outflow and continuum sources are 
given in Table 4. Disks are expected to be aligned orthogonal to the outflow directions, given that 


the protostar and disk together are thought to be responsible for outflow generation (Frank et al. 


2014). We have adopted a somewhat loose limit for the relative position angles for several reasons. 


First, outflow position angles are generally defined by eye and asymmetries in the outflow emission 
can lead to incorrect position angles. Second, outflows can wander with time (e.g., L1157-mm; 
Gueth et al.||1996) and have bends and kinks (e.g., L1448C; Hirano et al.||2010), adding uncertainty 


to the outflow axis. Lastly, the continuum emission from the protostellar envelope is entangled with 
that of the disk, potentially adding systematic uncertainty to the position angle of the continuum. 


3.1.1. Disk-like Structures 


Two sources in our sample show strong evidence of disk-like structures in the dust continuum: 
L1448 IRS2 (Figures [2] and [4]) and Per-emb-14 (Figure [2]). They are well-resolved in the dust 
continuum and have deconvolved position angles within 30° of being orthogonal to the outflow. 
The outflow directions of these sources are listed in Table 1 and known from both Spitzer IRAC 
imaging and/or CO outflow mapping. 


The disk-like structure toward L1448 IRS2 has a radius of ~300 AU and is at an angle of 71° 
with respect to the outflow in the B and C-array combined image (Figure [I]). The SMA Extended 
and Compact combined image (Figure [4]) also shows this structure. The B-array data alone do not 
appear symmetric and the emission is only extended toward the southwest at an angle of 68° with 
respect to the outflow. This extension does not appear to trace a discrete source and the SMA 
Extended image at 850 /un also exhibits this structure (Figure [4]). 


Per-emb-14, also known as NGC 1333 IRAS 4C (Smith et al. 2000), appears rather featureless 
in the C and B-array combined image (Figure [Ij), but the B-array data alone reveal a compact 
extended structure at an angle of 70° with respect to the outflow direction (Figure [2]). The radius 
of this structure is ~100 AU and resembles the expected structure of an edge-on protostellar disk. 


An outflow has not been detected toward this source in CO (Curtis et al. 2010 Plunkett et al. 


2013). However, a fan-shaped scattered light nebula is obvious in Spitzer maps as well as H 2 


emission knots extending to the east of the source position from which the outflow position angle 


°CASA stands for Common Astronomy Software Application (McMullin et al. 20071 and can be obtained from 
http://casa.nrao.edu 


















is derived; the IRAC imaging is shown in the Appendix. An edge-on orientation could also explain 
the lack of a CO outflow detection because the projected velocities would be low and likely confused 
with the cloud emission. 

In addition to Per-emb-14 and L1448IRS2, the sources L1448C, IRAS 03282+3035, and L1448 
IRS3C have deconvolved position angles that are within 30° of orthogonal to the outflow. Therefore, 
these sources could also be disk candidates, but they are not as well-resolved as Per-emb-14 and 
L1448IRS2; they will be discussed further in the following sections. 


3.1.2. Candidate Multiple Sources 

Our high-resolution data also reveal several candidate multiple sources with more than one 
continuum peak observed. Toward L1448 IRS3B we observe extended structure perpendicular to 
the outflow in the combined B and C-array images as well as the presence of another continuum 
peak. The B-array imaging alone reveals that there are at least two continuum sources present. The 
secondary source is located to the northwest by ~0 / . / 9. The secondary also appears to be marginally 
resolved, with a separation of ~0 / . / 2, slightly less than the minor axis of the beam ((//23). The 
nature of the secondary is not immediately clear since it is almost exactly along the axis of the 
outflow, but it may be a binary companion. The deconvolved position angle in Table 4 for the 
low-resolution image of L1448 IRS3B does not reflect the emission that is extended orthogonal to 
the outflow because the fitting routine attempts to fit both the main protostar and the companion 
along the outflow direction. 

IRAS 03292+3039 also has structured continuum emission that is extended in the north-south 
direction at an angle of 49° with respect to the outflow. There are no strong peaks standing out 
from this structure, even when imaging with higher resolution, but there are three weak sub-peaks 
separated from each other by ~0 / /4. 

3.1.3. Smooth Continuum Structures 

The sources IC348-MMS, L1448C, L1451-MMS have rather smooth continuum structures. 
They all have resolved envelopes extended over several resolution elements in the B and C combined 
imaging as shown in Figure [l] The higher resolution imaging from the B-array data alone does not 
reveal strong evidence of sub-structure and they are all consistent with smooth source structure at 
this resolution (Figure [2]). However, the deconvolved position angle of L1448C has an angle of 83° 
with respect to the outflow, perhaps indicative of a disk-like structure at the limit of our resolution. 
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3.1.4■ Asymmetric Resolved Structures 

Several sources show evidence of resolved structure at high-resolution, but these are not clearly 
disk-like in appearance and are often only extended toward one side. L1448 IRS3C (also called 
L1448 NW (Terebey & Pad gett||T 997}) appears resolved at an angle of 99° with respect to outflow, 
but only toward the north-east (Figure [2]). IRAS 03282+3035 is also extended at a slight angle 
(109°) with respect to the outflow direction, but only toward the north (Figure [2]). The Class I 
source L1448 IRS3A shows a roughly symmetric structure in the northwest and southeast directions 
(Figure [3]). The outflow direction from this source is uncertain (see Appendix and 
2006}, so this could be a disk-like structure, but it is uncertain without clear knowledge of the 
outflow direction. 

The strongest asymmetric structure is found toward IRAS 03292+3039 (Figures [I] and [ 2 ]). This 
source was also identified as a candidate multiple, due to its multiple peaks. However, its extended 
source structure stands out relative the other sources in the sample. We will discuss the nature of 
this source further in Section 5.2. 


Kwon et al. 


3.2. Known Multiple Systems 


Several systems in our sample were previously known to have wide components separated by 
more than l". L1448C is known to have a companion denoted L1448C-S by Jprgensen et al. (2006), 


also detected by Tobin et al. (2007) and Chen et al. (2013). IC348-MMS and IRAS 03282+3035 


were found to have wide binaries detected at millimeter wavelengths by Chen et al. (2013). Finally, 


L1448 IRS3 is comprised of three components (IRS 3B, IRS 3A, IRS 3C or NW; Looney et al 


2000 ). 


We have confidently detected all components of the L1448 IRS3 system, and we also find that 
L1448 IRS3B has a closer companion separated by only 1//9. We detect the Class I companion 
toward L1448C with a flux density of ~ 7 mJy, separated by 871 (~1900AU). However, we do not 
detect the companion toward IRAS 03282+3035, reported by Chen et al. (2013) with a separation 
of l'/5. We should have detected this source with a signal-to-noise ratio of 7.6. The non-detection 
of the companion may indicate that it is not real and it could be a feature extended along the 
outflow cavity. 

The wide binary (~15 w ) known as IC348-MMS2 is not detected because it lies near the edge 
of the CARMA primary beam, but it does have further detections at both submillimeter and far- 
infrared wavelengths ( |Palau et al. 2014). In addition, Rodriguez et al. (2014) detected IC348-MMS2 
and another source at 2 cm and 3.3 cm located ~ 3 7/ southwest of IC348-MMS, denoted JVLA 3a 
and IC348-MMS is JVLA 3b. JVLA 3a is undetected in our data and if it were emitting as a 
typical protostar, it should produce significant millimeter-wave emission; JVLA 3a source is also 


brighter than IC348-MMS at cm wavelengths. The Figures shown by Palau et al. (2014) indicate 
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that JVLA 3a may be coincident with the origin of the outflow emission at 4.5 ^m as well as the 
Spitzer 24 //rri source, while IC348-MMS (JVLA 3b) appears slightly offset. 


3.3. Molecular Line Observations 

We attempted to observe molecular line emission in the course of all the continuum obser¬ 
vations, but in nearly all cases the 13 CO and C 18 0 emission was not strongly detected due to 
observing 3 sources per track. Therefore, we obtained follow-up observations of the strong disk 
candidates L1448 IRS2 and Per-emb-14 with longer integration time in order to characterize their 
molecular line kinematics in the 13 CO and C 18 0 (J = 2 —> 1) transitions. The sources in L1448 
IRS2, L1448C, L1448 IRS2, L1451-MMS, and the three sources in L1448 IRS3 did have detections 
of their outflows in 12 CO emission (in addition to 13 CO for L1448 IRS2) and we discuss these 
observations in the Appendix. 


3.3.1. L1448 IRS2 

The presence of the large continuum structure toward L1448 IRS2 is suggestive of a disk, but 
molecular line observations are necessary to characterize its kinematics and determine its nature. 
The C 18 0 (J = 2 —> 1) integrated intensity maps and position-velocity diagrams are shown in 
Figure [5j The lower resolution image (top panels) is a combination of SMA Compact and CARMA 
C-array observations, showing a velocity gradient on ~5 W (1200 AU) scales in both the blueshifted 
and redshifted integrated intensity and position-velocity plot. The position angle between the peak 
blueshifted and redshifted emission is 70° east of north, making a 68° angle with respect to the 
outflow. Therefore, this velocity gradient could be interpreted as rotation. 

The SMA Compact and Extended map is shown in the bottom panel of Figure [5] with the 
respective position-velocity plot. The addition of the sensitive Extended array data shows that the 
velocity gradient direction changes to a position angle of 100° east of north and is now more in the 
direction of the outflow on ~2 W scales. There are several possibilities that could produce a velocity 
gradient along the outflow that we will discuss further in Section 5.2. Therefore, we cannot confirm 
rotation on scales less than ~300 AU toward L1448 IRS2. 

3.3.2. Per-emb-14 

We have also obtained higher sensitivity molecular line imaging toward Per-emb-14 (NGC 1333 
IRAS4C) in the C 18 0 (J = 2 —> 1) transition with CARMA in C-array; 13 CO (J = 2 —> 1) was 
not detected. The C 18 0 emission is weak, but we have made integrated intensity maps of the red 
and blue-shifted emission toward the source that are shown in Figure [6j The blue and red-shifted 
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emission are offset from each other, with the velocity gradient being perpendicular to the outflow 
and consistent with the axis of the resolved continuum source. The offset is only ~0 / . / 5, which is 
about 40% of the beam width; however, the velocity centroid accuracy is approximately equivalent 
to the beam width divided by the signal-to-noise ratio (Reid et al. 1988), 0'.'25 in this case. While 
the gradient is most likely real, the signal-to-noise is not high enough to determine if it is indicative 
of rotational support. 


3.4. Visibility Amplitude Profiles 

The well-sampled uv-plane provided by CARMA produces deconvolved images that reliably 
recover complex structure; however, the visibility amplitude data themselves still provide details 
of source structure that are not immediately apparent in the deconvolved images. The circularly 
averaged visibility amplitudes are plotted against projected uv-distance for each source and are 
shown in Figure [7j with uv-coverage out to ~600 kA. These data show a variety of structures from 
roughly constant amplitude out to 600 kA (0 / /33) to very little emission on baselines longer than 
~200 kA (scales smaller than 1"). 

L1448C and L1451-MMS exhibit the flattest visibility amplitudes. L1448C shows a linear 
decline in the log-log plots; the linearity is evident in L1448C out to the longest baselines. L1451- 
MMS shows the least amount of visibility amplitude decline of all the observed sources. 

Per-emb-14, IRAS 03282+3035, IC348-MMS, and IRAS 03292+3039 all decrease slowly until 
a scale specific to each source and then drop quickly. The scales are consistent with the size of 
the resolved structures apparent in the deconvolved images. L1448 IRS2 also had ~4 W diameter 
flattened structure, but the visibility amplitudes decrease more slowly than IRAS 03292+3039 
likely due to the structure being smaller in one dimension. L1448 IRS2 also shows variations at 
uv-distances > 100 kA, likely due to the resolved structure seen in the images on >1" scales. 

L1448 IRS3B showed clear evidence of multiplicity in its deconvolved image (Figure [2]). The 
visibility amplitudes flatten at >100 kA and then begin to drop at ~250 kA ((//8), the approximate 
separation of the candidate companion source. The visibility amplitude structures of L1448 IRS3C 
and L1448 IRS3A are difficult to interpret due to their proximity to L1448 IRS3B; moreover, both 
sources show abrupt drops in visibility amplitude at ~200 kA and ~ 400 kA. The visibility data 
will be further interpreted in conjunction with radiative transfer models in Section 4. 


3.5. Dust Continuum Masses 

The dust emission observed toward the protostars is directly proportional to the mass of the 
emitting dust, provided that the dust opacity, the temperature of the emitting dust, and the optical 
depth are known. At millimeter wavelengths the dust opacities are low enough such that optically 
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thin emission is a reasonable assumption. However, the dust opacity is only known to within a 


factor of 3 to 5 from models (e.g., Ossenkopf & Henning 1994 Ormel et al. 2011) and empirical 


estimates (e.g., Beckwith et al. 1990). Moreover, there can be opacity variations from source to 


source (Kwon et al. 2009 Chiang et al. 2012) and radial opacity variations in envelopes (Kwon 


et al. 2009) and disks (Perez et al. 2012). Since these sources all have protostars, there will also be 
a temperature gradient in the envelope. However for simplicity, we will assume a constant, average 
temperature of 30 K, appropriate for the environs within a few x 100 AU from the protostar. 

The mass of the emitting dust can be calculated with the relationship 

M = KxB x (T dust y (1) 

where D is the distance, k is the dust opacity at the observed wavelength, B(T dust ) is the Planck 
function, and T dust is the assumed temperature of the emitting material. We are adopting the dust 
opacities of Ossenkopf & Henning (1994) for this calculation (fti. 3 mm = 0.899 cm 2 g _1 , dust-only 


opacity), which are appropriate for the conditions of protostellar clouds. Finally, the total mass 
can be calculated by assuming a dust to gas mass ratio, the standard value of which is 1:100. The 
systematic uncertainty in the derived masses can be factors of several due to uncertainties in the 
dust opacity and dust temperature. 


We will use this method to estimate the mass of the emitting material from the observed flux 
density at 1.3 mm, assuming that we are only detecting thermal dust emission. To measure the flux 
densities, we will employ two methods. First, we will measure the flux density in the deconvolved 
images at both low and high resolution and use these flux measurements to calculate the mass. 
Second, we will measure the flux densities directly from the visibility amplitudes on 50 kA and 100 
kA scales. On these scales the envelope contribution is small and the residual contribution can be 
estimated from single-dish flux density measurements. 


3.5.1. Masses from Deconvolved Images 

To measure the continuum flux densities directly from the deconvolved images, we sum the 
flux above the 3a level within a 2.5" aperture around the protostar position. Thus, we only include 
the flux density bounded by our 3a image contours shown in Figures [l] and [2} We perform these 
measurements on both the lower resolution images and the high resolution images. The integrated 
flux densities, peak flux densities, RMS noise, and inferred masses are given in Table 4. 

The lowest mass measured in the low-resolution imaging is 0.017 Mq toward L1451-MMS 
and the largest mass is 0.48 Mq toward L1448 IRS3B. These mass measurements would include 
contributions from both the disk and envelope. The median and average low resolution masses are 
0.1 Mq and 0.13 Mq, respectively, with a standard deviation of 0.14 Mq. The masses from the high 
resolution data, on the other hand, have filtered more of the large-scale emission and may better 
probe the mass from the compact inner envelope and/or disk. The largest high-resolution mass is 
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0.19 Mq toward IRAS 03292+3039 and the lowest is 0.011 Mq toward L1451-MMS. The median 
and average high resolution masses are 0.067 Mq and 0.064 Mq, respectively, with a standard 
deviation of 0.054 Mq. 

The masses at both high and low-resolutions have variations over an order of magnitude, 
though the standard deviation of the high-resolution masses is over a factor of two less than the 
standard deviation at low-resolution. Moreover, there are no trends in mass (or flux density) with 
other protostar properties such as bolometric temperature (Tb 0 z) or bolometric luminosity (L boi). 


3.5.2. Masses from Visibility Amplitudes 


We also measured the masses directly from the visibility amplitude data that are shown in 
Figure [7J This enables us to determine better the spatial scales of the emitting material. We 
measured the flux densities at 50 kA and 100 kA, corresponding to ~4" (920 AU) and ~2" (420 
AU), respectively, and any emission must be coming from this or smaller scales. The flux densities 
and masses are given in Table 5. 

The median masses derived from the flux densities at 50 kA and 100 kA are 0.09 Mq and 0.07 
Mq respectively, with respective average masses of 0.1 Mq and 0.05 M 0 . The standard deviations 
of the visibility amplitude masses are lower than those measured from the deconvolved images, 0.08 
Mq and 0.03 Mq in units for 50 kA and 100 kA, respectively. 

Despite the isolation of emission from small-spatial scales, these masses still include any enve¬ 


lope flux that originates on these scales. Jprgensen et al. (2009) developed a method to correct for 
the residual envelope flux density at 50 kA. They determined that an envelope with a p oc R -1 " 5 
density profile will contribute at most 4% of its total emission at 50 kA; we verified this and also 
determined that at 100 kA the envelope contribution is only 2% of its total flux density. If the 
density profile were steeper, p oc R~ LS for instance, Jprgensen et al. (2009) found that the envelope 
could then contribute up to 8% to the flux density at 50 kA. 


The total envelope flux density for the sample in Jprgensen et al. (2009) was determined from 


SCUBA 850 pm data. For our sources, there were 1.1 mm Bolocam data taken toward all protostars 
in the Perseus cloud. We scaled the 1.1 mm flux densities to 1.3 mm, assuming (3 = 1.78 from 


Ossenkopf k Henning (1994); the flux densities at 50 kA are typically 10% to 25% of the Bolocam 


flux densities. We then subtracted the estimated envelope contribution of 4% and 2% of the scaled 
Bolocam flux density for 50 kA and 100 kA scales, respectively. 


As an independent check, we applied this method to the resolved disk of L1527 IRS (Tobin 


et al. 2013). The 870 gm and 3.4 mm visibility amplitudes at 50 kA were 314 mJy and 19.7 mJy 


respectively. The single-dish flux densities of L1527 at 870 pm and 3.3 mm are 1690 mJy and 33.3 
mJy respectively (Shirley et al. 2011). Assuming the single-dish flux densities only contribute 0.04 
of the 50 kA flux, the corrected 50 kA flux densities are 246 mJy and 18.4 mJy at 850 pm. and 3.4 



















14 - 


mm respectively. These values are close to the integrated flux densities from the highest resolution 


images of the resolved disk in Tobin et al. (2013). 


The median envelope-corrected masses at 50 kA and 100 kA are 0.052 Mq and 0.046 Mq 
respectively, with respective average masses of 0.068 Mq and 0.040 Mq. The standard deviations 
of the masses are 0.08 Mq and 0.03 Mq at 50 kA and 100 kA, respectively, with roughly the same 
amount of scatter as in Jprgensen et al. (2009). A histogram of disk masses at 50 kA and 100 kA 
is shown in Figure [8j We show the corrected and uncorrected masses at 50 kA and 100 kA versus 
T bo i in Figure i Our sample shows no sign of a correlation between the calculated masses and 


Tboh while the results of Jprgensen et al. (2009) may have shown a weak correlation. However, the 


T boi values we use are calculated with the inclusion of longer wavelength Herschel data (Sadavoy 
et al. 2014) and may be more accurate. For instance, Tboi = 47 K for L1448C in Sadavoy et al 


(2014), but Jprgensen et al. (2009) lists T boi = 77 K for this source and Enoch et al. (2009) gives 


T boi = 69 K. We also only cover a narrow range of T boi since we are focused on Class 0 sources, 


while Jprgensen et al. (2009) had equal numbers of Class 0 and Class I sources. 


4. Model Comparison 


The interpretation of the 1.3 mm continuum emission can be aided by comparison to radiative 
transfer models of axisymmetric envelopes, disks, and envelopes with embedded disks. We have run 


a small grid of models using the Hyperion code (Robitaille 2011). This is a Monte Carlo radiative 


transfer code that calculates radiative equilibrium and generates high signal-to-noise images of 
dust emission using ray tracing. With these models, we aim to determine if the data are consistent 
with envelope-only models or if a compact structure is required. We ran three varieties of models: 
disk-only models, rotating-collapsing envelopes (CMU envelopes; Cassen Sz Moosman||1981 Ulrich 


1976), and power-law envelopes; both types of envelope models are run with and without disk 


components. 

The models are run with three variations of disk mass 0.0 Mq, 0.01 Mq, and 0.1 Mq and 
with disk radii of 10 AU, 30 AU, 50 AU, 100 AU, 300 AU, and 500 AU. The disks have a surface 
density profile oc R _1 and a scale height of 10 AU at a radius of 100 AU; the surface density profile 
and scale height are chosen to be consistent with a viscous accretion disk in vertical hydrostatic 
equilibrium (Shakura Sz Syunyaev 1973). The envelopes all have masses of 5.25 Mq, a 1 Lq central 
star, an outer radius of 10,000 AU, and an outflow cavity with a half-opening angle of 10°. The 
envelope mass corresponds to a CMU model with a mass infall rate of 5 x 10~ 5 Mq yr -1 and 
this was chosen to be representative of a young protostellar system dominated by the mass of the 
envelope. The inner radius of the envelope and disk are equivalent to the dust destruction radius, 
adopted to be where the dust temperature is 1600 K. For the CMU envelopes, we assumed that the 
disk radius was equivalent to the centrifugal radius (Re); the radius where material is rotationally 
supported. The power-law envelopes have either a p oc R -15 or R~ 2 spherical density profile; these 


density profiles are representative of free-fall collapse and the singular isothermal sphere (Shu 1977) 
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(or alternatively the Larson-Penston solution (Larson 1969)), respectively. 


The models were set up to generate output images at 1.3 mm with 5 AU pixels and we 
Fourier transformed the images using the MIRIAD task fft. We azimuthally averaged the Fourier 
transformed images to compare with the observed visibility amplitudes. The normalized visibility 
amplitudes from the envelope-only models (0.0 M 0 disks) and the disk-only models are shown in 
Figure [lOj As expected, the CMU models with increasing Re have more rapidly declining visibility 
amplitudes and a power-law envelope with p oc R -2 has a more slowly decreasing visibility amplitude 
than the R -15 envelope and CMU envelopes. The disk-only models have flat amplitudes for disks 
with small radii, but disks with radii > 100 AU show a factor of 5 decline at uv distances > 100 kA 
and a 100 AU disk has a factor of 5 decline at 500 kA. The disk-only models also show variation 
in visibility amplitudes due to the disks with smaller radii being optically thick, resulting in less 
emergent flux. 


Disks embedded within envelopes are shown in Figures 0 and [T2j The addition of the disk 

components does add some flattening to the visibility amplitude profiles, but they are always 

declining. The model with Re = 500 AU does, however, have a small upward bump. The disks 
also change the ratio of the visibility amplitudes at small uv-distances relative to large uv-distances 
(large scales to small scales, respectively). For assumed envelope masses of 5.25 M 0 , a 0.01 M 0 

disk does not dramatically affect the visibility amplitude profiles, but a 0.1 M 0 disk has a more 

profound effect. We also notice that a 100 AU, 0.1 M 0 disk has a flatter visibility amplitude profile 
between 100 kA and 400 kA than disks with radii between 10 AU and 50 AU. This is due to the 
small, massive disks being optically thick as mentioned in the previous paragraph. 


A qualitative comparison of the model visibility amplitude profiles to the data reveal that 
many sources have profiles that decline more slowly than the models of envelopes with embedded 
disks. Furthermore, some sources look more like the disk-only models than disk+envelope models 
(e.g., Per-emb-14, IRAS 03282+3035, IC348-MMS). 

A more quantitative comparison of the models and observations is shown in Figure [13J where 
we plot the ratio of visibility amplitude at 25 kA to the visibility amplitude at 250 kA. This plot 
gives a sense of the relative contribution of the envelope and compact component (modeled as a 
disk) on spatial scales that are an order of magnitude different, 10" (~2300 AU) and 1" (~230 
AU) scales. These scales are chosen for both physical and practical reasons; 25 kA is the minimum 
scale that we could measure accurate visibility amplitudes due to our chosen array configurations. 
Furthermore, a 10" (~2300 AU) scale is likely to be directly associated with gravitational collapse 
and formation of the protostar, while a disk is expected to form on scales of about 1" (~230 AU) 
or smaller. Lastly, at uv-distances greater than 250 kA, not all sources had sufficient signal-to-noise 
to accurately measure their visibility amplitude. 

The disk-only models are shown as plus signs in Figure [l3| with an increasing ratio as the disk 
radius increases and the power-law envelopes without disks are plotted as the large symbols. We 
also show the models with increasing disk radii within a given envelope (triangles and squares). An 
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envelope with a 10 AU, 0.1 Mq embedded disk has large ratios because the disk is optically thick, 
and the ratios decrease with disk radius until R = 100 AU and then the ratios begin to increase 
again as the disk itself becomes resolved. These tracks with disk radii can be raised or lowered by 
changing the ratio of disk mass to envelope mass. 

The key result shown in Figure [l3| is that for a 25 kA to 250 kA ratios less than ~8 (the ratio of 
the p oc R~ 2 envelope, the large diamond in Figure [TTi| a compact component of some size and mass 
is required to explain the visibility amplitudes. The most likely form of such a compact component 
is a protostellar disk. This criterion is fulfilled for 6 out of 9 Class 0 protostars: L1451-MMS, 


also demonstrates that massive (~ 0.1 Mq), embedded disks with radii < 30 AU may not be 
apparent toward protostars at 1.3 mm because the 25kA to 250 kA visibility amplitude ratio are 
not significantly lower than the envelope-only case. This happens because emission from the small, 
massive disks being optically thick masking much of the emitting material. 


Per-emb-14, L1448 IRS2, L1448C, IC348-MMS, and IRAS 03282+3035. Furthermore, Figure 13 


5. Discussion 

The results presented here are currently the largest collection high-resolution 1.3 mm contin¬ 
uum observations of Class 0 protostars in a single star forming region. Moreover, our data all have 
comparable resolution and sensitivity, a feature lacking from previous studies. 


5.1. Disk Formation 


The formation of disks during the protostellar phase has been thought to occur as a con¬ 
sequence of angular momentum conservation during protostellar collapse (|Ulrich 1976; Cassen &; 


Moosman 

1981 

Terebey et al. 

1984) 


the protostellar phase, dependent entirely on the angular momentum inherited from the infalling 
envelope. Hydrodynamic simulations, without magnetic fields, readily form large disks during col¬ 


lapse (e.g., Yorke & Bodenheimer 1999). These disks are massive enough to be gravitationally 
unstable, forming spiral arms and fragments (e.g., Boley|[2009 Vorobyov||2010 Kratter et al.pOlO 


Zhu et al. 2012) 


However, magnetic fields, which had been shown to potentially slow rotation during the pre¬ 


collapse phase (e.g., Basu 1998), were shown by Allen et al. (2003) to also significantly slow rotation 
during collapse and suppress the formation of a rotationally-supported disk (Galli et al. 2006). 


These studies were verified by Mellon & Li (2008) and Hennebelle & Fromang (2008); disk formation 


was prevented even in the presence of very weak magnetic fields. More recent studies showed that 
by including non-ideal MHD effects in simulations, initially small rotationally-supported disks could 
form (Dapp <fc Basu 2010). The infalling material would still undergo magnetic braking, but the 


magnetic fields would be dissipated in the high-density material close to the protostar, enabling 
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disk formation of a small rotationally supported disk. The size of the rotationally-supported region 


will grow with time, but is expected to remain < 10 AU throughout the Class 0 phase (Dapp et al 


2012). Machida et al. (2008) also simulated the formation of disks with non-ideal MHD and found 


that massive 100 AU disks could form while in the Class 0 phase. A study by Joos et al. (2012) 


explored less ideal initial conditions, with misaligned magnetic fields and rotation axes. These 
simulations showed that even in the ideal MHD limit disks could form if the magnetic fields are 
misaligned with respect to the rotation axis. 

The misaligned magnetic field scenario may be plausible because the TADPOL survey of 


magnetic field morphologies (Hull et al. 2013, 2014) found no systematic alignment of magnetic 
fields and outflow axes (presumed to reflect the rotation axis). Moreover, two of the four known 


Class 0 systems with rotationally-supported disks (L1527 IRS and VLA 1623; Tobin et al. 2012 


Murillo et al. 2013) were found to have magnetic fields perpendicular to the outflow axis on ~1000 


AU scales and down to the scale of the disk in L1527 IRS (Segura-Cox et al. 2015). The sources 


L1448C, L1448 IRS3B, and L1448 IRS2 had observations in the TADPOL survey. L1448C only had 
a few detected magnetic field vectors, but they are at a ~45° angle to the outflow, L1448 IRS3B 
had vectors that are at an angle of ~90°, and L1448 IRS2 had vectors that are both aligned and 
misaligned with respect to the outflow. 

To summarize, theory predicts two different scenarios for disk formation, that depend on the 
included physics and the initial conditions, specifically the treatment and importance of magnetic 
fields. However, recent work with misaligned magnetic fields and non-ideal MHD (e.g., Tomida 


et al. 2015) seems to have reduced the significance of the ‘magnetic braking catastrophe.’ 


5.2. Evidence for Class 0 Disks 


Thus far, rotationally-supported disks have been detected toward four Class 0 protostars: 
L1527, VLA 1623, RCrA IRS7B, and HH212 MMS ( Tobin et al.||2012 Murillo et al.||20l3 Codella 


et al. 2014 Lindberg et al. 2014). L1527, VLA 1623, and RCrA IRS7B all have both extended 


dust emission perpendicular to the outflow direction and Keplerian rotation signature, while HH212 
MMS is only resolved kinematically. It is unclear if the small number of observed Class 0 disks 
reflects a true paucity due to the physics of disk formation or if it is simply a lack of observations 
with high enough resolution and sensitivity to detect Class 0 protostellar disks (or disk candidates). 
Moreover, given that the emission from the disks is entangled with that of the envelope, they are 
more difficult to directly resolve than disks around Class I and Class II sources with much less (or 
no) envelope emission. 


To increase the likelihood of detecting disks embedded within infalling envelopes, our sample 
is slightly biased toward sources that might have orientations that are within 30° of edge-on and 
well-defined outflow axes from either infrared scattered light or CO emission. This is because edge- 
on disks should stand out better against the surrounding envelope due to the higher column density 
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through the disk midplane. While our data are very sensitive and have a resolution of ~0 / . / 3 (70 
AU), this is barely at the limit where we would expect to detect evidence for disks around Class 0 
protostars. Indeed, for a source two-beams across, the diameter would have to be at least ~140 AU. 
The spatial resolution achieved by our survey is a factor of ~ 1.6 more coarse than the observations 
toward L1527 IRS in Tobin et ah (2012). If L1527 IRS was at the distance to Perseus, it would 
only be marginally resolved in our observations. 


5.2.1. Sources with Evidence for R > 100 AU Disks 


Despite the difficulties posed by limited spatial resolution, we do detect two strong candidate 
disks toward L1448 IRS2 and Per-emb-14 (NGC 1333 IRAS4C). The dust continuum structure 
around L1448 IRS2 is very large, ~300 AU in radius. Per-emb-14 is more compact and has an 
apparent radius of ~100 AU. 

The new and archival kinematic data for L1448 IRS2 did not enable a clear detection of 
Keplerian rotation around the protostar. Near the protostar, the velocity gradient in C 18 0 only 
differs from the outflow position angle by 30°. The dust continuum emission at higher resolution is 
extended toward the southwest and Tobin et al. (in prep.) detects a binary source at the end of this 


dust extension. Takakuwa et al. (2014) showed that infall from a circumbinary disk, shepherded by 


the binary sources could produce a radial velocity gradient along the outflow direction. Thus, the 
rotation signature on >200 AU scales could be that of a circumbinary disk and the shift in velocity 
gradient direction could reflect material flowing through the disk toward the binary sources. On 
the other hand, a velocity gradient along the outflow direction is suspicious given that the outflow 


can entrain the ambient envelope material (Arce & Sargent 2006). Thus, we merely suggest radial 


transport of material through the disk as a possibility to produce the velocity gradient along the 
outflow. 

Per-emb-14 also appears to be a promising disk candidate from its resolved dust continuum 
structure. Our C 18 0 data detected evidence of a velocity gradient in the expected direction for 
rotation. However, higher sensitivity will be required to verify this velocity gradient and determine 
if this rotation is Keplerian. 

In addition to these two most promising R > 100 AU disk candidates, IRAS 03292+3039 
has one of the brightest, largest, and most puzzling structures. While the position angle of this 
source is not close to being orthogonal to the outflow, the combined imaging from C and B- 
configurations shows a ~2.5 W (575 AU) diameter structure. Then at the highest resolution, the 
emission is rather constant across the source with a brightness temperature of ~8K, indicating 
that the dust may be close to being optically thick, in which case the brightness temperature 
would reflect the temperature of the emitting material. In addition to the size of this structure, 
strong evidence of rotation on ~3 /7 scales was shown byiSchnee et al. (2012) in C 18 0 (J = 2 —> 1) 
observations. Furthermore, the velocities of the outflow emission traced by CO (J = 2 —> 1) 
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indicate that this source is viewed at an intermediate inclination, possibly 40° (Yen et al. 2015). 
The C 18 0 (J = 2 —> 1) data was further analyzed by Yen et al. (2015), finding that a disk as large 
as 850 AU could form in this system. 

The large resolved structure, intermediate inclination, and rotation on scales comparable to 
the continuum size are evidence that this source may also be a good disk candidate, even though the 
continuum emission is not extended orthogonal to the outflow direction. The mass of the continuum 
source is also large, ~0.2 M 0 at the highest resolution and the protostar mass is estimated to be 


~0.3 Mq (Yen et al. 2015). If the continuum source is indeed a rotationally-supported, the large 
ratio of disk to protostar mass could mean that the disk would be gravitationally unstable. However, 
future observations with higher resolution will be needed to verify if the continuum structure is a 
rotationally supported disk as well as a more robust kinematic measurement of the protostar mass. 


5.2.2. Sources with Evidence for R < 100 AU Disks 


While there is evidence for large disk-like structures around a few protostars, the results are 
not so clear for the remaining sources. The sources L1448C, IRAS 03282+3035 and L1448 IRS3C, 
have deconvolved position angles that are 7°, 19°, and 9° from being orthogonal to the outflow, 
respectively. However, the continuum emission is not obviously disk-like in the images, but only 
higher-resolution/sensitivity data will be able to determine this for certain, in addition to molecular 
line observations. 


The analysis of visibility amplitude ratios enables us to shed further light on the evidence for 
disks (or at least compact dust structures) on scales smaller than 100 AU. Figure 13 shows that 
sources with 25 kA to 250 kA amplitude ratios less than ~8 require a contribution from a compact 
component, possibly a disk. A ratio of 8 is the value for a spherical envelope with a radial density 
profile p oc R~ 2 and no compact density structure, the steepest we might expect in idealized star 
formation models; the ratio of an R -15 density profile is ~30. The ratio expected for a 500 AU 
disk itself is ~5. Thus, the ratio between 0 and 8 will depend on both the radius of the disk and 
the amount of mass in the disk relative to the envelope. 


L1451-MMS, Per-emb-14, L1448C, IC348 MMS, L1448 IRS2, and IRAS 03282+3035 all have 
a 25 kA to 250 kA ratio less than 8 (6 out of 9 Class 0s in the sample). L1448 IRS3C was discussed 
previously as having resolved emission nearly orthogonal to the outflow, but it does not have a ratio 
less than 8. This could result from its close proximity to L1448 IRS 3B, given that its visibility 
amplitudes have unexplained dips at 200 kA and 40 kA. IRAS 03292+3039 (discussed in the previous 
section) also has a ratio > 8, but this could be caused by the emission from the disk and inner 
envelope being optically thick out to > 1" scales, leading to the emission being resolved-out. 


L1448 IRS3B appears to have extended envelope emission orthogonal to the outflow, encom¬ 
passing the two sources. The detected emission is not significantly larger in angular extent than 
that of L1448 IRS2 and IRAS 03292+3039, but the emission is not centered on either of the con- 
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tinuum sources. Kinematic observations of L1448 IRS3B by Yen et al. (2015) show that there is a 
velocity gradient consistent with rotation on 3" scales, but higher resolution kinematic observations 
will be needed to determine the nature of this structure. 


5.2.3. Outlook 

The results from this small survey have shown that at least two (Per-emb-14 and L1448 IRS2), 
possibly three (including IRAS 03292+3039) show evidence for resolved disks with radii > 100 AU. 
The remaining sources have evidence for compact structure within their envelopes that could be 
evidence for a disk, but higher resolution continuum and kinematic observations are necessary to 
characterize them. Despite the small sample, it is encouraging for studies of Class 0 disks that 3 of 9 
Class 0 sources have evidence for resolved disk-like structure and many of the others have compact 
emission extended perpendicular to their outflows or visibility amplitude ratios that suggest a disk- 
like component. While we cannot prove that the disk candidates are rotationally supported, it is 
important to highlight that we are detecting an abundance of structure on scales less than 500 AU 
toward Class 0 protostars, possibly a signature of protostellar disks. 


5.3. Class 0 Disk Masses 


Assuming that we are probing emission from the protostellar disks, we have been able to 
calculate their masses to compare with measurements of other Class 0 sources and more evolved 
Class I and II objects. The median masses of the compact components are ~0.05 Mq (values 
corrected for estimated envelope emission) for measurements taken at both 50 kA and 100 kA. This 
mass is an order of magnitude larger than the characteristic mass of Class II disks (~0.005 Mq, 
Andrews &; Williams|2005 ) and 5x larger than the median Class I disk mass measured in Jprgensen 


et al. (2009). The Class 0 disk masses measured in Jprgensen et al. (2009) did have some sources 


with even larger masses than we find, which may support a scenario of typically higher disk masses 


in the Class 0 phase. A caveat of this comparison is that while 

Jprgensen et al. (2009) assumed 

the same dust opacities as our study, Andrews & Williams 

(2005 

) assumed = 0.035 (dust 

plus gas opacity), about a factor of two larger than the 

Ossenkopf &; Henning (1994|) model at 850 


fin i (a factor of two larger opacity decreases the calculated mass by the same factor). However, a 
larger opacity might be more appropriate for Class II disks, but more importantly the gas to dust 
mass ratio in Class II disks might be significantly less than 100 (Williams Sz Best 2014). The gas 
to dust ratio uncertainty is potentially more significant than the opacity uncertainty in the case of 
Class II disks, while the gas and dust may still be well-mixed in the Class 0 phase. 


The main uncertainty in masses derived in this analysis is the true envelope density structure. 
Numerous studies have shown that Class 0 envelopes may have a variety of density profiles; the 
most typical profiles are steeper than p oc R _1 5 


(Looney et al. 2003, Chiang et al. 


Kwon et al. 
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2009, Chiang et al. 2012; Tobin et al. 2015). This is important, because steeper envelope density 


profiles result in greater contributions of envelope emission at 50 kA and 100 kA, leading to over 
estimated disk masses. 


5.4. Multiplicity 


Our sample is not large enough to make strong revisions to the current multiplicity results 


toward Class 0 protostars, which was recently found to be 0.64 by Chen et al. (2013). We have 


discovered one new close companion (L1448 IRS3B separated by ~210 AU) and we believe that the 
300 AU companion toward IRAS 03282+3035 is spurious. Thus, the multiplicity fraction as a whole 
does not change. However, we did sensitively sample smaller spatial scales than previous studies 
and only confidently find 1 new clear candidate companion (L1448 IRS 3B). IRAS 03292+3035 
may be another candidate, but its multiplicity is uncertain. 


Maury et al. (2010) had suggested a lack of multiplicity on scales between 150 AU and 430 AU 
(corrected for the updated distance to Perseus), but our new binary detection is separated by ~210 


AU and the sample of Chen et al. (2013) survey was quite incomplete at those scales. A larger, 


less-biased sample will better illuminate the multiplicity frequency as a function of separation; such 
a project is being carried out now with the VLA (Tobin et al. 2015 in prep.). 


6. Summary 

We have presented a CARMA 1.3 mm survey of 9 Class 0 protostars and 2 Class I protostars 
in the Perseus molecular cloud. This is one of the largest high-resolution (~0 / /3, 70 AU) samples of 
1.3 mm data taken toward Class 0 protostars (thus far). We also include kinematic follow-up for 
two sources (Per-emb-14 and L1448 IRS2). The main results from the survey can be summarized 
as follows. 

1) We detect three strong Class 0 disk candidates toward L1448 IRS2, Per-emb-14, and IRAS 
03292+3039. These systems are not yet kinematically confirmed to be Keplerian disks. L1448 
IRS2 shows evidence of rotation on scales > 200 AU, with a change in velocity gradient direction 
on scales < 200 AU. Per-emb-14 has an indication of rotation on the scale of the detected disk¬ 
like structure, but the spectral line data have low signal-to-noise and we cannot determine if the 
rotation is Keplerian. IRAS 03292+3039 does show evidence for rotation on 500 AU scales, but 
smaller scale measurements are not yet available. 

2) A variety of resolved structures are detected within our continuum data, aside from the 
disk candidates listed in 1), ranging from apparent massive inner envelopes or disks (L1448 IRS 3B 
and IRAS 03292+3039) to marginally resolved structures on < 200 AU scales (IRAS 03282+3035, 
L1448 IRS3C). Finally, three sources do not appear to have resolved structure on scales < 200 AU 
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(L1448C, IC348 MMS, L1451-MMS). Moreover, 5 out of 9 sources in the sample have deconvolved 
position angles that are within 30° of orthogonal to the outflow, possibly an indication of compact 
structure consistent with a disk. 


3) Comparison of the observed visibility amplitude ratios to radiative transfer models of disks 
and envelopes enables us to infer that at least 6 of 9 sources in our sample require a compact 
component, possibly represented by a disk (L1451-MMS, Per-emb-14, L1448C, IC348 MMS, L1448 
IRS2 and IRAS 03282+3035). The radiative transfer models also show that emission from embedded 
disks with masses ~0.1 Mq and radii < 30 AU will be optically thick and have lower levels of 
emission at 1.3 mm relative to larger disks with the same mass. 


4) A candidate companion to L1448 IRS3B is detected with a separation of Cf/9 (~210 AU) 
and the companion itself may be comprised of two sources separated by ~0 / /2. IRAS 03292+3039 
may have multiple components, but the multiple peaks detected do not significantly stand out from 
the resolved structure. We do not detect a companion source toward IRAS 03282+3035 that had 
been detected by Chen et al. (2013). 


We conclude that there is evidence for the formation of large disks in the Class 0 phase, but 
Class 0 disks likely have a range of radii and masses that depend on the structure, kinematics, 
and possibly magnetic field properties of their parent cores. Sub-arcsecond resolution imaging is 
crucial to characterizing the structure of Class 0 disks and observations with ALMA and the VLA 
at wavelengths longer than 1.3 mm will likely be necessary to characterize emission from small, 
massive disks. 
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A. Outflow Maps 

Our observations included the 12 C0 and 13 C0 (J = 2 —> 1) transitions that, when detected, 
enabled us to examine the outflow structure toward some of our protostellar targets at resolutions 
that are higher than typically achieved for outflow emission. 


A.l. L1448 IRS3 

The L1448 IRS3 region is of particular interest because of the complex, overlapping outflows 


in the region (Kwon et al. 2006) and there have not been higher angular resolution observations 


toward these sources. We show the outflow channel maps in Figure 14 from a combination of 


our higher-resolution data and the CO data observed in CARMA D-array by Hull et al. (2014), 


gaining better sensitivity to extended structures. The sources are marked with crosses and we can 
clearly identify the outflow of L1448 IRS3B, the southern most source. There appears to be a 
contribution from the L1448C outflow in the southern parts of the map. We also can identify, for 
the first time, an outflow associated with the northern-most source, L1448 IRS3C. For this source, 
we mainly detect the redshifted side of the outflow and the blueshifted side is observed in only a few 
velocity channels. Nevertheless, both the blue and red-shifted sides of the outflow have a cone-like 
morphology that can be traced back to an origin at the position of L1448 IRS3C. The asymmetric 
brightness of the blue and red-shifted outflow components may result from L1448 IRS3C being 
located near the edge of the L1448 cloud and having less entrained material on the blue-shifted 
side. 


Kwon et al. (2006) had associated what we now believe to be the red-shifted side of the L1448 


IRS3C outflow with L1448 IRS3A, the more evolved Class I protostar in the region. However, the 
bipolar pattern of the CO emission seems to clearly originate at the position of IRS3C. We do see 
some redshifted CO emission near L1448 IRS3A that seems to be extended toward the northeast; 
however, we cannot clearly identify a position angle for this outflow. Scattered light morphology 


shown by Tobin et al. (2007) may suggest that the outflow has a slight north-west position angle, 
but the complexities of the emission from the three protostars within 30" of each other make such 
an identification difficult. 















-24- 


A.2. L1451-MMS 


We detected CO (J = 2 —> 1) toward L1451-MMS, the candidate first hydrostatic core identi¬ 
fied by Pineda et al. (2011), who also found a compact CO outflow in their lower resolution data. 
We confirm the presence of the CO outflow with our higher-resolution data as shown in the channel 
maps in Figure 15 Moreover, we also find that the continuum emission is slightly offset from the 
line connecting the blue and red-shifted outflow; this slight offset was also seen in the data from 


Pineda et al. (2011). The offset may be related to the outflow being driven into an inhomogeneous 


medium. This offset is also in the direction of the 3cr extension seen in the continuum emission. 
The detection of the outflow with our high-resolution data demonstrates the compactness of the 
outflow, otherwise it would have been resolved-out. 


A.3. L1448 IRS2 


We detect the outflow toward L1448 IRS2 in CO (J = 2 —> 1), as shown in Figure 16, with 
excellent sensitivity. The high resolution and sensitivity is achieved with a combination of SMA 
Compact and Extended configurations as well as CARMA C and B configurations. The outflow 
is clearly offset from the 1.3 mm continuum peak by about 1" and the outflow cavities appear 
asymmetric, opening wider to the southwest than the north east. This asymmetry is evident in 
both the blue and red-shifted lobes of the outflow. The asymmetry in the opening angle could reflect 
the structure of the ambient medium in the direct vicinity of the protostars, with the envelope being 
less dense to the southwest and enabling the cavities to open more widely. Furthermore, the outflow 
would appear more symmetric if the driving source was shifted southwest a few arcseconds, along 
the extended dust continuum emission. 


The combined sensitivity of the CARMA and SMA maps also enabled 13 CO (J = 2 —> 1) to be 
resolved in the outflow toward L1448 IRS2 (Figure [l7|. The 13 CO emission traces more compact 
structure located near the continuum source, but the asymmetry and spatial offset evident in 
the 12 CO emission are still present. Moreover, the location of the blue-shifted intensity peak is 
coincident with the blue-shifted intensity peak in C 18 0, evidence that the compact C 18 0 may be 
affected by the outflow. 


A.4. L1448C 


Lastly, we were able to detect the compact inner outflow toward L1448C in CO (J = 2 —> 1) 
with only the CARMA B and C-array data. We were able to obtain a good detection with such 
little time on source due to L1448C having extremely bright CO emission. Maps with superior 


sensitivity and comparable resolution do exist (e.g., Hirano et al. 2010), but we include the map 
for the sake of completeness. 
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A.5. Per-emb-14 


Per-emb-14 does not currently have a CO outflow detection, possibly due to the source being 
edge-on. However, the Spitzer IRAC images of the source reveal a scattered light nebula at the 
position of the source, along with shock features extended along the inferred blueshifted side of the 
outflow. These data were taken as part of Spitzer GO program 30516 and were reduced with the 


same methods described in Tobin et al. (2010b) and are more sensitive than the data presented 


by Gutermuth et al. (2008). We infer that the eastern side of the outflow is blueshifted because 


it is brighter and it is known that the geometrical effects will cause the redshifted side to be more 
extincted and have a fainter nebula (Whitney et aL|[2003). 
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Fig. 1.— Continuum images at 1.3 mm of the full sample, emphasizing the structure of the sources 
on the scale of a few arcseconds. Images are ordered from top to bottom by their amount of resolved 
structure; the most resolved structures are at the top. The contours in all images are [-3, 3, 6, 9, 
12, 15, 20, 25, 30, 35, 40, 45, 50, 60, ...] x a, where a is denoted in each panel and in Table 4. 
The beam size is shown in the lower right corner of each image and exact beam sizes are given in 
Table 4. The blue and red arrows denote the direction of the blueshifted and redshifted outflows, 
respectively. 
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Fig. 2.— Same as Figure 1, but with emphasis on sub-arcsecond structures. All plots are zoomed-in 
relative to Figure 1, except IRAS 03292+3039. 
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Fig. 3.— Continuum images at 1.3 mm of the neighboring Class I sources that were within the 
primary beam: L1448 IRS3A (top panels) and L1448C-S (bottom panels). The contours in all 
images are [-3, 3, 6, 9, 12, 15, 20, 25, 30, 35, 40, 45, 50, 60, ...] x a, where a is denoted in each 
panel and in Table 4. The beam size is shown in the lower right corner of each image and exact 
dimensions are given in Table 4. The blue and red arrows denote the direction of the blueshifted 
and redshifted outflows, respectively. 











- 32 - 



RA offset (arcsec) 


RA offset (arcsec) 


Fig. 4.— L1448 IRS2 images at 1.3 mm (Extended & Compact) and 850 /im (Extended-only) from 
the SMA. The 1.3 mm image traces similar structure as the CARMA 1.3 mm images in Figures 
1 and 2. The 850 /.fin image at higher resolution still traces larger-scale structure, but does show 
the asymmetry evident in the higher resolution 1.3 mm image. The contours in the 1.3 mm images 
start at 7a, increasing by 2a until 15<r at which point they increase by 5er; the 850 //rri contours 
increase in the same sense, except that they start at 3a. At 1.3 mm and 850 /im a = 1.2 mJy 
and 3.5 mJy, respectively. The beam sizes are 1.11" x 0.92" and 0 / . / 65 x (y/54 at 1.3 mm and 850 
gm, respectively. The blue and red arrows denote the direction of the blueshifted and redshifted 
outflows, respectively. 
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Fig. 5.— Integrated intensity images of red and blue-shifted C 18 0 (J = 2 —> 1) emission overlaid 
on 1.3 mm dust continuum imaging and position-velocity diagrams taken from the region enclosed 
by the parallel dashed lines (right). The top panels are from CARMA C and B-array and SMA 
C 18 0 Compact data, with C//5 tapering applied; the lower resolution data show a clear velocity 
gradient on 5 W (~1200 AU scales) at an angle of 70° with respect to the outflow position angle. 
The bottom panels show the higher resolution SMA Compact and Extended array observations of 
C 18 0; the direction of the velocity gradient changes to have an angle of only 30° with respect to 
the outflow position angle on scales less than 2" (~500 AU). The contours in the CARMA and 
SMA imaging (top) are [-6, 6, 9, ...]x a where sigma is 0.46 K and 0.48 K for the red and blue 
shifted emission respectively. The contours for the SMA-only imaging are [-5, 9, 13, 17] x a where 
sigma is 0.68 K and 0.62 K for the red and blue shifted emission respectively. The larger beams 
are for the line data and are T/82 x T/48 and T/31 x T/29 for the low and high resolution images 
respectively; the smaller beams are for the continuum data, the same as in Figures 1 and [4] The 
blue and red arrows denote the direction of the blueshifted and redshifted outflows, respectively. 
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Fig. 6.— Per-emb-14 (NGC 1333 IRAS4C) with CARMA C 18 0 overlaid on the B and C-array 
dust continuum image. The integrated intensity of the red and blue-shifted emission is suggestive 
of rotation on the scale of the disk candidate, but the signal-to-noise is far too low for us to 
determine if the rotation is consistent with Keplerian. The contours are ±3 a and 5cr, where a = 
1.74 K and 1.35 K for the red and blue-shifted maps respectively. The larger beam drawn in the 
lower right corresponds to the C 18 0 data and is T/19 x l 7 /11; the smaller beam is for the continuum 
data and is the same as in Figure 2. The blue and red arrows denote the direction of the blueshifted 
and redshifted outflows, respectively. 
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LI448 IRS3B 



LI 448 IRS3C 



L1448C LI448 IRS2 




Fig. 7.— Binned visibility amplitude plots versus projected uv-distance at 1.3 mm. The visibility 
data are averaged within 25 kA bins. The sources show a variety of structures from flat visibility 
amplitudes toward L1451-MMS, power-law decline in L1448C, high scatter for L1448 IRS2 indicat¬ 
ing resolved structure, steep declines at a variety of uv-distances (Per-emb-14, IRAS 03292+3039, 
IRAS 03282+3035), and evidence of multiplicity L1448 IRS3B. The dotted line is the expected 
visibility amplitude from noise alone. 
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50 kX Masses 



Log Moss (M 0 ) 
100 kX Masses 



Log Mass (M©) 


Fig. 8.— Histogram plot of disk masses inferred from the flux densities at 50 kA (top) and 100 
kA (bottom). Masses that have been corrected for estimated envelope contribution (solid line) and 
with no correction (dashed line) are shown. The median disk masses corrected for envelope emission 
are 0.052 Mq and 0.046 Mq at 50 kA flux and 100 kA, respectively. The respective median disk 
masses that are not corrected for envelope emission are 0.09 Mq and 0.07 Mq, respectively. 
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100 kX Mass, Corrected 
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Fig. 9.— Plot of (corrected) 50 kA disk masses versus L b 0 i (top left) and T bd (bottom left); the 
same plots but for 100 kA masses are plotted on the right. There is no apparent trend in the inferred 
disk masses as a function of L m or T boi- Jprgensen et al. (2009) also found inconclusive evidence 
for Class 0 disk masses to trend with these parameters. The lack of correlation indicates that the 
disk properties (at least mass) may not depend specifically on evolutionary state and luminosity, 
but rather initial conditions. 
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CMU Envelope, 0.0 M sun disk, Variable R c Spherical Envelope p a R p , No disk 




No Envelope, 0.01 M sun disk 



No Envelope, 0.1 M su „ disk 



Fig. 10.— Visibility amplitudes for envelope models with and without disk components. The upper 
left shows a standard rotating collapse model (CMU) without a disk component. The upper right 
shows power-law envelopes with density profiles proportional to R~ 15 and R” 2 0 with no embedded 
disk. The lower left shows disk-only models with sizes between 10 AU and 500 AU with a total 
mass of 0.01 M@. The lower right shows disk-only models having a total mass of 0.1 Mg. The 
visibility amplitudes of the disk-only models are not normalized in order to show that for fixed 
mass, the disks with small radii have high optical depths thereby reducing the emergent 1.3 mm 
flux. 
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CMU Envelope, 0.01 M sun disk, Variable R c 



CMU Envelope, 0.1 M slin disk, Variable R c 



Fig. 11.— Visibility amplitudes for CMU Envelope models with disk components. The left panel 
shows CMU envelope models with an embedded disk of sizes between 10 AU and 500 AU with a 
total mass of 0.01 Mq. The right panel shows CMU envelope models with disks having a total 
mass of 0.1 Mq. The envelopes have a total mass of 5.25 Mq. 
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Spherical Envelope p a R 1 ' , 0.01 M sun disk 




Spherical Envelope pa R '' , 0.1 M sun disk 




Fig. 12.— Visibility amplitudes for Power-law Envelope models with disk components. The upper 
left panel shows p oc R ~ 15 envelope models with embedded disks of sizes between 10 AU and 500 
AU with a total disk mass of 0.01 Mq. The upper right panel is the same as the upper left, except 
that the disks have a total mass of 0.1 Mq. The lower left panel shows p oc R~ 2 0 envelope models 
with embedded disks having a total mass of 0.01 Mq. The lower right panel is the same as the left, 
except that the disks have a total mass of 0.1 Mq. The envelopes have a total mass of 5.25 Mq. 
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0-0 M sun disk 



Fig. 13.— Visibility amplitude ratio plots for the visibility amplitudes at 25 kA and 250 kA; 
corresponding to 10" (2300 AU) and 1" (230 AU) scales, respectively. The sources are plotted as 
the solid circles and the models are plotted as the crosses (disk only), triangles (CMU models), 
squares (p oc R~ 15 ), and diamonds (p oc R~ 20 ). This hgure (a) shows models without a disk 
component, but still a rotationally-flattened region for the CMU case. Figure 13b shows models 
with a 0.01 Mq disk component, and Figure 13c shows models with a 0.1 M@ disk component. These 
plots imply that sources with a 25 kA to 250 kA ratio less than ~8 require a disk or another compact 
component, or a density profile steeper than p oc R~ 2 '°. Envelope models with an embedded disk 
can be made to have lower ratios by decreasing the envelope mass, emphasizing the disk component. 
The centrifugal/disk radii are labeled at 10 AU, 100 AU, and 500 AU for the CMU and disk models 
such that the variation of the visibility amplitude ratios can be traced. The F^(25kA) values for 
the disk-only models have been scaled by the factors noted in each plot such that the points for all 
radii are visible. 






F„(25 kX)/F,(250 kX) 


100 


0.01 



a CMU 
□ p = — 1.5 

o p = —2.0 
+ Disk only x 


tn 


m 


o 

t n 
< 
cr 

FQ 

500 AU£ 
+ - 

'+ ° 
T o 


5.0 


CN 

CO 

Y 

00 



\f 

\|_ 100 AU 

. J 10 AU. + + 

1 [DnresoTvfecr 



J. 


100 

F.(25 


Fig. 13b. 


LI448 IRS 3C 


M sun disk 



kX) (mjy) 


000 







sun 


disk 



(mjy) 


000 









-46- 



3 h 25 m 37?557?0 36?5 36?0 35?5 35^25 m 37f537?0 36?5 36?0 35?5 33^05 m 37f557fO 36?5 36?0 35?5 35^05 m 37f537fO 36?5 36?0 35?5 35?0 

RA (J2000) RA (J2000) RA (J2000) RA (J2000) 


Fig. 14.— Outflows in the L1448 IRS3 region. The sources are marked with crosses; the northern 
most source is L1448 IRS3C, the middle source is L1448 IRS3A, and the southern-most source is 
L1448 IRS 3B. The blue-shifted emission is shown in Figure 14a, while the red-shifted emission is 
shown in Figure 14b. The blue-shifted outflow from IRS 3B becomes apparent at -4 km s” 1 and 
the blue-shifted outflow from IRS3C becomes visible at 1 km s _1 . The emission visible south of 
IRS3B, may be from L1448C. The red-shifted emission is apparent both IRS 3C and IRS3B at 7 
km s -1 out to 17 km s^ 1 . There is some red-shifted emission associated with IRS3A between 7 
km s -1 and 9 km s _1 , but it is unclear if this is outflow emission from this source. The contours 
are [-5, 5, 10, 15, 20, ...] X a and cr = 0.1 Jy beam^ 1 ; the beam is 2 / . / l x l'!7. The blue and red 
arrows denote the direction of the blueshifted and redshifted outflows, respectively. 





































Dec (J2000) w Dec (J2000) w Dec (J2000) w Dec (J2000) 


- 47 - 



3 h 25 m 37?537?0 36?5 36?0 35?5 35?0 

RA (J2000) 


Fig. 14b. 
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Fig. 15.— Compact outflow from L1451-MMS. The outflow is detected to be in the same direction 
as observed by Pineda et al. (2011) and with a similar offset to the east of the 1.3 mm continuum 
peak. The contours are ±[2, 3, 4, 5] x a and <7=0.1 Jy bearn^ 1 ; the beam is Cf/92 x 0 , . , 72. The 
blue and red arrows denote the direction of the blueshifted and redshifted outflows, respectively. 
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Fig. 16.— Integrated intensity maps of 12 C0 emission toward L1448 IRS2 from the combined imag¬ 
ing of CARMA and SMA datasets, overlaid on the CARMA 1.3 mm continuum map (grayscale). 
The emission maps reveal a rather asymmetric appearance of the outflow cavities, opening wider to 
the southwest in comparison to the northeast. Furthermore, the outflow vectors, which are centered 
on the peak of the continuum emission, seem as if they would be better centered slightly to the 
southwest, along the extended continuum emission. The contours are [-6, 6, 9, 12, 15, 18,...] x 
a and er=2.47 K, 2.55 K for the red and blue-shifted integrated intensity maps, respectively. The 
beam for the CO emission is 1"42 x l" 13 and the beam for the continuum is (//91 x R/68. The 
blue and red arrows denote the direction of the blueshifted and redshifted outflows, respectively. 
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RA (J2000) 


Fig. 17.— Integrated intensity maps of 13 CO emission toward L1448 IRS2 from the combined imag¬ 
ing of CARMA and SMA datasets, overlaid on the CARMA 1.3 mm continuum map (grayscale). 
In the case of L1448 IRS2, the 13 CO emission is tracing the dense, lower velocity outflow. The 
positional offset of the blue-shifted emission is consistent with that of the C 18 0 emission shown in 
Figure [5j The contours are [-3, 3, 6, 9, 12, 15, 18,...] x a and <7=1.14 K, 1.44 K for the red and 
blue-shifted integrated intensity maps, respectively. The beam for the CO emission is \".‘2 x T/06 
and the beam for the continuum is (//91 x (y/68. The blue and red arrows denote the direction of 
the blueshifted and redshifted outflows, respectively. 
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RA (J2000) 


Fig. 18.— Integrated intensity maps of 12 C0 emission toward L1448C from the CARMA B and 
C-array data, overlaid on the CARMA 1.3 mm continuum map (grayscale). The contours are [-4, 4 
6, 9, 12, 15, 18,...] x a and a= 13.922 K, 15.39 K for the red and blue-shifted integrated intensity 
maps, respectively. The beam for the CO emission is (//90 x C//75 and the beam for the continuum 
is Cf/58 x Cf/41. The blue and red arrows denote the direction of the blueshifted and redshifted 
outflows, respectively. 
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Right Ascension (J2000) 


Fig. 19.— Spitzer IRAC 3.6 pm, 4.5 pm, and 8.0 pm imaging of the NGC 1333 IRAS4 region. The 
blue and red vectors mark the outflow position angle adopted for Per-emb-14 (NGC 1333 IRAS4C). 
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Table 2. CARMA 1.3 mm Observation Log 
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Sources Config. Date Track Length Antennas Central Frequencies T sys Gain Calibrator/Flux Density Flux Calibrator 

(UT) (hrs) (GHz) (K) (Name, Jy) 
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Table 4. 1.3 mm Continuum Properties 
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and position angles are determined using the imfit task in CASA. The error in the relative angle between the deconvolved source position angle and outflow position angle assumes no error in 
the outflow position angle. 



Table 5. 1.3 mm Continuum at 100 kA and 50 kA 
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